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We show that, Landau level mixing in two-dimensional quantum dot wave functions can be 
taken into account very effectively by multiplying the exact lowest Landau level wave functions 
by a Jastrow factor which is optimized by variance minimization. The comparison between exact 
diagonalization and fixed phase diffusion Monte Carlo results suggests that the phase of the many- 
body wave functions are not affected much by Landau level mixing. We apply these wave functions 
to study the transition from the maximum density droplet state (incipient integer quantum Hall 
state with angular momentum L — N(N — l)/2) to lower density droplet states (L > N(N — f)/2). 



I. INTRODUCTION 

Advances in miniaturization and processing techniques 
of semiconductor systems have made it possible to create 
nanostructures containing only a small number of elec- 
trons IHEl IE El- These quantum dots, also called ar- 
tificial atoms, provide a crucial testing ground for many 
important quantum effects in low-dimensional systems. 
The system parameters can be tuned by controlling the 
dot geometry, electrostatic gate voltage, and by applying 
a magnetic field. 

Experimental techniques such as capacitance spec- 
troscopy 0, Q or gated transport spectroscopy 0, Q al- 
low experimentalists to precisely control the number of 
electrons in a two-dimensional parabolic confining poten- 
tial, and map their electronic properties as a function of 
the magnetic field. Such systems offer very rich quantum 
many-body effects, providing a link between two differ- 
ent branches of physics: In the limit of zero magnetic 
field the system acts like an atom, i.e., it exhibits shell 
structure and obeys Hund's first rule. On the other hand, 
for relatively modest magnetic fields the system exhibits 
incipient quantum Hall physics. A large number of tran- 
sitions induced by the magnetic field can be observed 
experimentally. The transitions can be understood from 
Hartree-Fock theory or from density functional theory (at 
low magnetic fields even the simple constant interaction 
model J'} suffices) but at high magnetic fields an accurate 
treatment of many-body effects is important. 

The aim of this work is to study in detail the tran- 
sitions between many-body states for moderately large 
magnetic field values where the angular momentum of 
the system increases beyond that of the so called max- 
imum density droplet (MDD) regime Q], corresponding 
to the filling factor v = 1 in quantum Hall physics. Most 
of the previous work on the transition between MDD to 
lower density states (LDD) was performed in the lowest 
Landau level (LLL) approximation and /or was restricted 
to spin polarized electrons |8ll9l ll0lllllll2lll3| . The LDD 
states are not always predicted to be fully spin polarized, 
both in the absence of Landau level (LL) mixing Q and 



in the presence of LL mixing |14| . 

Here, we combine exact diagonalization and quantum 
Monte Carlo (QMC) techniques [H [H E3 to obtain 
a very accurate description of many-body states beyond 
the MDD. Linear combinations of determinants obtained 
by exact diagonalization within the LLL approximation 
serve as the starting point for building trial wave func- 
tions for our variational Monte Carlo (VMC) and dif- 
fusion Monte Carlo (DMC) calculations. The method 
described here is very systematic, can be applied to any 
eigenstate, and, is more efficient than exact diagonaliza- 
tion with LL mixing. Trial wave functions obtained by 
the exact diagonalization method have the right sym- 
metry properties, including conservation of total angular 
momentum L, total spin S, and separability of center of 
mass and relative motion, which allows us to study center 
of mass excitations 18]. We apply our method to inves- 
tigate the many-body states involved in the MDD-LDD 
transition. We find that for the N — 4 electron system, 
the pair densities of the ground and excited states close 
to the MDD-LDD transition have a well defined square 
symmetry, while those for N = 6 states have little struc- 
ture. We find that the MDD-LDD transition acts as a 
convergence point for special values of (L, S) which cor- 
respond to the magic numbers . The redistribution of 
electronic charge is rather smooth in the absence of the 
Zeeman effect, but, when a Zeeman term is included in 
the Hamiltonian it becomes more abrupt. 



II. MODEL AND METHOD 

A. Hamiltonian 

We consider an N— electron two-dimensional quantum 
dot with a circularly symmetric parabolic confinement 
potential. Within the effective mass approximation, and 
neglecting finite thickness effects, the Hamiltonian in an 
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external magnetic field, B, is 




where to* is the effective mass of electron, A is the 
vector potential, luq is the parameter characterizing the 
parabolic confinement potential, e is the dielectric con- 
stant, g* is the effective (/-factor, [iq is the Bohr mag- 
neton and s z j is the z-component of the spin of the j th 
electron. 

We use standard materials constants corresponding to 
GaAs, to* — 0.067mo, dielectric constant e = 12.4 and 
g* = —0.44. Results will be given in effective atomic 
units (Ti = e 2 /e = c = to* = 1): the effective Bohr radius 
is (Zq = 9.79373 nm, and the effective Hartree is H* = 
11.8572 meV. We choose the confinement parameter hu>o 
to be 0.28 H* or 3.32 meV, which is within the typical 
experimental range [3- 0- HI • 

B. Single-particle states 

Ignoring the interaction term, the singleparticle Fock- 
Darwin energies of this Hamiltonian are |2(| 

e nrl = (2n r + \l\ + l)hu - (2) 

where lu c = eB/m*c is the cyclotron frequency, u) — 
+ n r is the radial quantum number and I 

is the angular quantum number (n r , \l\ = 0,1,..., oo). 
Equivalently, the energy spectrum can be expressed as 
a sum of two harmonic-oscillators [tj , 

e nm = (n + l/2)huj + + {m + l/2)htu^, (3) 

where u>± = u) ± lu c /2, n is the Landau level index, 
n,m — 0, 1, oo and I — —n, —n + 1, oo. The re- 
lationship between these two sets of quantum numbers 
is I = m — ri, n r = min(n, to) or n = n r + 

n , + max(0, —I), m = n r + = n r + max(0, 1). 

At zero B the energies depend only on the principal 
quantum number, 2n r + \l\ = n + to, whereas at infi- 
nite B they depend only on the Landau level index n. 
The corresponding eigenfunctions \nma) (<r is the spin 
index) are used as the basis set in our exact diagonaliza- 
tion and QMC calculations of interacting dots. At zero 
magnetic field the system has degeneracies and therefore 
a shell structure, as in the case of real atoms. When a 
magnetic field is applied, the degeneracy is broken and 
several transitions occur due to single-electron level cross- 
ings as described by Eq. [3 Eventually all the electrons 
are in lowest Landau level orbitals and each of these or- 
bitals is doubly occupied. In a noninteracting electron 



picture no more transitions are predicted as the mag- 
netic field is increased further, but, experimentally they 
are observed @ . These can be understood from Hartree- 
Fock or density functional theory. The intra Landau level 
splittings get smaller with increasing B and the exchange 
interaction favors parallel spins, leading to spin flips, un- 
til eventually all the spins are parallel and the N lowest 
angular momentum states are occupied singly. This is 
the maximum density droplet (MDD) state. As B is in- 
creased further, Hartree- Fock [2lJ and density functional 
calculations [22, |23| predict an edge reconstruction due 
to some electrons jumping from low angular momentum 
states to higher angular momentum states because the 
Coulomb interaction overwhelms the parabolic external 
potential. These lower density droplet (LDD) states were 
observed experimentallv|24| but quantitative discrepan- 
cies with experiments exist. Since several determinants 
of LLL noninteracting orbitals are exactly degenerate it 
is expected that an accurate treatment of correlation is 
important. This is the range of B that we study in this 
paper. 

C. Many- particle states and Exact diagonalization 

The many-body wave functions may be chosen to be 
eigenstates of the total angular momentum L = L z , total 
spin S 2 and the z-component of total spin S z . In ad- 
dition, the dynamic symmetry of an isotropic parabolic 
potential leads to the separability of center of mass and 
relative motion and the existence of two center of mass 
operators [l8| . C+ and C- that commute with the H, L, 
S 2 and S z , and the above operators: 
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C- = 
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Determinants of Fock-Darwin orbitals arc eigenfunctions 
of L and S z . Thus, the many-body Hamiltonian in a de- 
terminantal basis is block diagonal, with the blocks cor- 
responding to different (L, S z ) values. Also, we note that 
LLL determinants are eigenfunctions of C+ with eigenval- 
ues C+ = 0. The size of the subspaces can be further re- 
duced by employing linear combinations of determinants, 
known as configuration state functions (CSFs), that are 
eigenstates of S 2 and C_ . The diagonalization of S 2 and 
C_ is done as explained in Ref.l)I8T). The matrix elements 
of the Coulomb interaction can be determined either by 
numerical integration or analytically [2?l l28l l2*flj . 

For a state with a given total angular momentum L, 
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and for a given number of Landau levels, n, only single- 
particle states with angular momentum l min < I < l max 
contribute, where Z m i n = —n. For fully spin polarized 
LLL wave functions, i max = L — (N — 1)(N — 2)/2 and 
for lower spin states it is higher than that. As the num- 
ber of Landau levels is increased, Z max increases, both 
because the remaining N — 1 electrons can lower their 
angular momenta by occupying more than one Landau 
level and because the n th Landau level has states with 
I > —n. Empirically we find that determinants contain- 
ing orbitals with / > 2L/N + 1 contribute little to the 
wave function and consequently we could limit ourselves 
to orbitals with I < 2L/N + 1. However, in most cases 
we avoid doing this because the wave functions are no 
longer exact eigenstates of the operator C- . 

D. Quantum Monte Carlo 

In this work, we use both the variational Monte Carlo 
(VMC) and the diffusion Monte Carlo (DMC) meth- 
ods [ia,llall3 to improve the exact diagonalization wave 
functions and energies obtained within the LLL. First, we 
introduce variational parameters by multiplying the LLL 
wave function by a Jastrow factor and optimize the Jas- 
trow parameters by minimizing the variance of the local 
energy ^30|. For a given state (L, S, S z , C_), the form of 
our trial wave functions is 

Ncsf 

' " ~ = JeeJedJeed ^ ' " ~J ( 6 ) 

i 

where J ee , J ec i, J ee d are the electron-electron, electron- 
dot, and electron-electron-dot Jastrow terms respec- 

L S S C 

tively. The CSF's ^ £ : z ' ~ are obtained by diagonaliz- 
ing the S 2 and C_ matrices. Multiplying by the Jastrow 
factor does not alter the L, S 2 and S z symmetries of the 
wave function, aside from introducing a completely neg- 
ligible amount of spin contamination [3l| (not an exact 
eigenstate of S 2 ), arising from the different cusp condi- 
tions for parallel- and antiparallel-spin electrons. The 
Jastrow factor reduces the statistical error in both VMC 
and DMC for a given number of Monte Carlo steps. It 
also lowers the VMC energy, but leave the fixed-phase 
DMC energy (described later) unchanged because the 
phase of the trial wave function is not altered. 

The linear coefficients ctj are obtained by diagonaliz- 
ing H. They can also serve as additional free parame- 
ters that can be reoptimized in VMC. However, we have 
found that, at sufficiently large magnetic fields, there is 
no gain in energy from reoptimizing them, and the exact 
LLL wave function multiplied by the optimized Jastrow 
factor provides an excellent trial wave function. This 
is in contrast to the situation in the absence of a mag- 
netic field, where the reoptimized linear coefficients in the 
presence of the Jastrow factor are considerably smaller 
in magnitude than the original ones from exact diago- 
nalization. This can be understood as follows. Quan- 



tum chemists distinguish between static correlation (also 
known as near-degeneracy correlation) and dynamic cor- 
relation. Within the quantum Monte Carlo literature 
it is well understood that static correlation is most ef- 
fectively incorporated into the wave function by using 
a linear combination of nearly degenerate determinants, 
whereas dynamic correlation is most efficiently treated by 
a flexible Jastrow factor. In the case of dots in sufficiently 
large magnetic fields the separation of energy scales (in- 
tra Landau level versus inter Landau level) is very large 
and the distinction between static and dynamic correla- 
tion becomes sharp. 

For any given state, the determinantal coefficients a^, 
for a LLL wave function, need not be recalculated when 
the system parameters m* , ujq, B and e are varied, be- 
cause the a.i are independent of these system parameters. 
This is because all the LLL determinants with a given L 
have the same expectation value for the single-particle 
Hamiltonian. So the single-particle terms in the Hamil- 
tonian contribute a multiple of the unit matrix and only 
the interaction part of the Hamiltonian need be diago- 
nalized. Hence, the LLL wave function is independent of 
e, and, depends on m*, luq and B only through a scaling 
of the linear dimensions by the factor yJTJmFuj. Conse- 
quently these Jastrow-LLL wave functions are not only 
a very accurate choice but also a very efficient choice 
for QMC calculations. Of course this does not apply to 
low angular momentum L states because it is impossible 
to construct them using only LLL orbitals. Note also 
that the independence of the oti on the system param- 
eters makes the Fock-Darwin orbitals the most efficient 
choice, as opposed to using say Hartree Fock or density 
functional orbitals. 

The DMC method improves upon the VMC energy by 
projecting onto the lowest energy wave function that has 
the same nodes (or the same phase in the case of complex 
wave functions) as the trial wave function. The resulting 
fixed-node or fixed-phase energies are upper bounds to 
the true energies and consequently the accuracy of the 
calculations depends on the quality of the nodes or the 
phase of the trial wave functions. We employ an efficient 
implementation of DMC described in Ref. |35|. Since we 
use complex trial wave functions, we employ the fixed 
phase approximation |32| , which has previously been ap- 
plied to quantum dot systems |33l l34| : starting from a 
variational wave function ^t, a new function / = 
is introduced, where 4> is the ground state of the iV-body 
Hamiltonian subject to the constraint that it have the 
same phase as V^t- The integral form of the Schrodinger 
equation in imaginary time (t — > it) then can be written 
in terms of / as: 

/(R,t + r) = J dR'G(R,R',T)f(R',t), (7) 

where R represents the coordinates of N electrons, 

G(R,R',t) = * T (R)(R|e- Tj& |R')/*T(R'), is a Green 
function that is explicitly known only in the short-time 
limit (t — ► 0). Hence, the Eq.([7|) can be solved by apply- 



4 



ing G repeadetly until the desired projection is obtained. 
Since the function / is used as a probability density for 
Monte Carlo sampling, it must be real and positive which 
is equivalent to the assumption that the phase of $ is the 
same as that of the trial wave function Wt- After equili- 
bration, the energy can be calculated as an average over 
the sampled points R m : 

1 M 

^-g^W. (8) 

m 

where El is the real part of the local energy, 

E L (R) = Re|* T (R)" 1 H* T (R)| . (9) 

We note, however, that for an external potential with 
circular symmetry it is also possible to use the fixed node 
approximation: In the subspace of the good quantum 
number L, the Hamiltonian can be written as 

H = y( ^-p* + \m*^ + - T — I - \^cL 
^ \ 2m* v ° 2 3 e f^rtj J 2 

(10) 

which is purely real. We can ignore the last term in H 
since it just gives a trivial shift in the energies of the 
various L states. For the remaining part of H, states 
with angular momentum L and — L are degenerate and 
so we can consider their sum, which is real. It is possi- 
ble to expand these real LLL wave functions in terms of 
real orbitals, but the number of determinants increases 
greatly. For instance, a MDD state can be written as a 
single determinant of complex orbitals or a sum of 2 _1 
determinants of real orbitals. So, it is more efficient to 
use complex wave functions within the fixed phase ap- 
proximation. 

III. RESULTS AND DISCUSSION 

We now apply the techniques described above to quan- 
tum dots similar to those studied experimentally 0,0, 
in magnetic fields ranging from 3 — 8 Tesla bringing the 
dots in the vicinity of MDD-LDD transition. Since the 
Zeeman term just results in a trivial shift of the differ- 
ent S z states relative to each other, we will ignore the 
Zeeman term unless explicitly stated otherwise. 

A. Comparison between Exact diagonalization and 

QMC 

In this subsection, we compare results obtained from 
exact diagonalization and QMC calculations with differ- 
ent levels of approximation. Fig. ^ shows the energies 
of N = 4 and N = 7 fully spin polarized electrons at 
B = 8 T and B = 5 T respectively. The solid line repre- 
sents the results obtained by exact diagonalization within 



the LLL approximation, from which we build QMC trial 
wave functions. The dotted lines are the exact diagonal- 
ization results with 2 LLs included. The DMC results, 
using the Jastrow-LLL wave functions, are represented 
with the crosses. The numbers represent C_ , with in- 
dicating that there is no center of mass excitation. As 
seen from Fig.^ DMC energies are very accurate, better 
than those from a 2 LL diagonalization. This is partic- 
ularly clear in Fig. IHb) where the magnetic field is only 
5 T since LL mixing is more important at lower mag- 
netic fields. For 4 spin-polarized electrons there is no real 
computational advantage to using QMC over exact diag- 
onalization, but, for spin unpolarized dots and/or higher 
numbers of electrons the QMC method becomes more 
efficient, because the number of determinants needed is 
dramatically smaller. For instance, in Fig.QJb) the QMC 
calculation for L = 29 using 21 determinants in the trial 
wave function clearly gives a better energy than a 2 LL 
diagonalization using 28637 determinants. These results 
also verify the appearance of magic numbers, i.e. partic- 
ular angular momenta, L* , at which a deep dip occurs, 
followed by a center of mass excitation at L* + 1 0, ^| . 
Also, we note that for N — 7, within the LLL approxi- 
mation exact diagonalization makes an incorrect predic- 
tion for the ground state which underscores the need to 
go beyond the LLL approximation near the MDD-LDD 
transition. 
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FIG. 1: Energy as a function of total angular momentum L, 
for (a) N = 4 at B = 8 T, (b) TV = 7 at B = 5 T, obtained 
by exact diagonalization in LLL approximation (solid line), 
exact diagonalization including 2LLs (dotted line), and DMC 
using full Jastrow term (cross marks). The numbers above 
the line represent the C_ quantum number. 
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TABLE I: Comparison of energies for different L values (N = 
4, S — 2, B = 8 T) obtained from exact diagonalization 
method within 1,2, and 3 LL approximations (ED), VMC 
including only the electron-electron part of the Jastrow factor 
(e-e VMC), and DMC using a fully optimized Jastrow factor 
(DMC). The energy differences between DMC and the 3-LL 
exact diagonalization, shown in the last column, are small, 
with the DMC energies being better in most cases. 



L 


ILL ED 


2LL ED 3LL ED 


p-p VMC 


DMC 


sYJ IV I V ; 1 J 1 J 


ft 
U 


^ ^0.99 




0.0400 




r. 





—29 x 10~ 4 


7 


tJ.uUUU 


5.4159 


5.4100 




5 


4071 


on v 1 n -4 


8 


5.4764 


5.3582 


5.3511 


5.3517 


5 


3491 


—20 x 10 -4 


9 


5.3609 


5.2688 


5.2641 


5.2634 


5 


2622 


-19 x 10~ 4 


10 


5.2544 


5.1733 


5.1705 


5.1738 


5 


1702 


-3 x 10~ 4 


11 


5.3182 


5.2370 


5.2343 


5.2376 


5 


2339 


-4 x 10" 4 


12 


5.2683 


5.2070 


5.2042 


5.2064 


5 


2041 


-1 x 10~ 4 


13 


5.2825 


5.2272 


5.2227 


5.2239 


5 


2218 


-9 x 10~ 4 


14 


5.1860 


5.1401 


5.1394 


5.1417 


5 


1394 


-0 x 10~ 4 


15 


5.2498 


5.2038 


5.2031 


5.2058 


5 


2033 


2 x 10 -4 


16 


5.2718 


5.2282 


5.2272 


5.2301 


5 


2275 


3 x 10~ 4 


17 


5.2965 


5.2563 


5.2548 


5.2574 


5 


2550 


2 x 10~ 4 


18 


5.2194 


5.1948 


5.1946 


5.1953 


5 


1946 


x 10 -4 


19 


5.2831 


5.2586 


5.2583 


5.2590 


5 


2584 


1 x 10 -4 



while the QMC results are not particularly superior to 
3 LLs results, agreement still remains excellent as seen 
in Table HU Since both VMC and DMC calculations are 
performed within the fixed phase approximation, we con- 
clude that exact LLL wave functions have very accurate 
phases. 
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The energies from 3 LL diagonalizations are indistin- 
guishable from the DMC energies on the scale of the 
plots, so we present them as tables. In Table Q] we show 
the energies for N = 4, S = 2, B = 8 T and in Table|TI]for 
N = 4, S = 0, B = 5 T. We have also included the results 
obtained by including only the electron-electron part of 
the Jastrow factor. The electron-electron Jastrow recov- 
ers most of the missing energy by keeping electrons apart, 
in agreement with the VMC calculations of spin polarized 
electrons in Ref. |3fJ For spin polarized electrons DMC 
gives better energies than 3 LL calculations especially 
for the lower L values. For higher L, electrons occupy 
larger orbits, making the electron-electron Jastrow less 
important. Thus, the effect of LL mixing is less impor- 
tant and exact diagonalization results agree with QMC 
results within a few standard deviations (« 1 x 10~ 4 
H*) of the QMC results. For spin unpolarized electrons, 



2 4 1 i i ' i ' i ' i ' i 

3 4 5 6 7 8 

B[Tesla] 

FIG. 2: Energy as a function of magnetic field for different 
many-body states for N = 4, (a) without Zeeman term, (b) 
with Zeeman term. Lines with symbols represent many-body 
states which are closely involved in the MDD-LDD transition, 
and they are identified by their (L,S) values in (a), and by 
their (L, S, S z ) values in (b). The locations of the maxima in 
the pair densities are indicated by solid dots, and, open dots 
represent the fixed electrons. 

Having established the accuracy of the form of the 
wave function, in this subsection, we study in detail the 
MDD-LDD transition for N = 4 and N = 6. Fig.0shows 
the energies as a function of magnetic field for different 
many-body states (L, S) (S z is also given in Fig. H[b)) 
closely involved in MDD-LDD transition. In Fig. |2Ja), 
the Zeeman splitting is not taken account, so the states 
are 2S + 1-fold degenerate. It should be noted that, in 



these calculations only the states with angular momen- 
tum N(N - l)/2 < L < N(N + 1)/2 are considered. For 
the fully spin polarized case the MDD-LDD transition 
involves only the L = N(N - l)/2 and L = N(N + l)/2 
states 0, 0| , but LLL exact diagonalization energies 
indicate that other angular momentum states between 
these values are expected to become ground states for 
spin unpolarized or partially polarized electrons Q . The 
excitation spectrum in Fig.[5{a) reveals a very interesting 
structure, involving several many-body states in a small 
region near the MDD-LDD transition which occurs at 
B rs 5.5 T: While these states (shown by solid lines with 
symbols) don't seem to have a regular predictable (L, S) 
pattern, they correspond to the magic number sequences 
of different S states for N = 4 [H SI E3, Strikingly, 
the MDD-LDD transition acts as a convergence point 
for these particular (L* , S*) values. On the other hand, 
when the Zeeman effect is taken into account (we con- 
sider an effective Lande factor g* = —0.44), higher spin 
states are favored and the separation between the solid 
lines with symbols and dashed lines become less clear. 

An interesting property of these states can be revealed 
by plotting their pair-density p(ri,r 2 ) (probability den- 
sity of two electrons being at ri and r^), with the ref- 
erence electron r 2 fixed at a density maximum. As 
schematically shown in Fig.|2Ia), the lower energy states 
all have the same structure, with electrons arranged in 
a square. In fact, square symmetry corresponds to the 
minimum energy electronic arrangement of classical point 
charges 0] . The higher energy states (dashed lines) have 
their electrons arranged in a centered triangle, except for 
the two center of mass excited states (L, S, C_) = (7, 2, 1) 
and (9, 0, 1). These two states show 3 and 5 peaks respec- 
tively in their pair-densities instead of the more intuitive 
4-peak structure (including the reference electron) one 
would expect for a N = 4 dot. We have also observed 
similar structures in other center of mass excitations for 
higher L values. 

In Fig. |3| we study the ground state pair densities at 
magnetic fields between 5 — 6 T as the system goes from 
the MDD state to the LDD states. Without Zeeman 
splitting, there is a smooth electronic redistribution in- 
volving the states (6,2), (8,0), (9,1), and (10,0). with 
increased localization of electrons. Note that, in this case, 
the state (10, 2) never becomes the ground state in agree- 
ment with the results of Ref . |l"8T) . When the Zeeman 
term is included, low spin states are suppressed and the 
system becomes fully polarized. As a result, the MDD- 
LDD transition is much more abrupt, going directly from 
the state (6,2) to the state (10,2). An abrupt charge 
redistribution in post-MDD transitions was observed ex- 
perimentally also . The dominant determinant in the 
(10, 2) state has LLL orbitals with I = 1,2, ...N and con- 
sequently the density at the center of the dot is small. 

In Fig. 0] we plot the energy spectrum for N = 6. 
As in the N = 4 case, we can distinguish between two 
sets of eigenstates when the Zeeman term is neglected 
(Fig.4(a)). We have only labeled the states that are 
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FIG. 3: (color online) Evolution of VMC ground state pair 
densities as a function of magnetic field for N = 4. (a) With- 
out Zeeman effect (b) Including Zeeman effect. 



closely involved in the MDD-LDD transition, marked 
by symbols. However, unlike the N = 4 case, we have 
found that, at the confinement strength considered here 
(uiq = 3.32 mcV), the pair densities of N = 6 dots near 
MDD-LDD transition have little structure (except for the 
(21, 3) state which has a small hexagonal-shaped ripple). 
In Fig. we plot the radial densities of the ground states 
that participate in the MDD-LDD transition, namely, 
(15,3), (18,0), (19,1), and (21,1). These ground state 
transitions are in agreement with those of Ref. [l4j except 
for the last transition where they observe a large region 
where the (20, 2) state occurs, followed by (21, 3). In our 
calculations, as seen in Fig. the (21, 3) state occurs 
only if a Zeeman term is included. As a result of the 
Zeeman splitting, the charge redistribution during the 
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FIG. 4: Energy as a function of magnetic field for different 
many-body states for N = 6, (a) without Zeeman term, (b) 
with Zeeman term. Lines with symbols represent many-body 
states which are closely involved in the MDD-LDD transition, 
and they are identified by their (L, S) values in (a), and by 
their (L,S,S Z ) values in (b). 



MDD-LDD transition becomes more abrupt, as in the 
N = 4 dot. 



IV. CONCLUSIONS 

To summarize, by combining exact diagonalization and 
QMC methods, we have obtained highly accurate many- 
body wave functions and energies for circular parabolic 
dots in the range of magnetic fields corresponding to 
the MDD-LDD transition. In this regime, the dot goes 
through several angular momentum and spin transitions, 
and one does not expect Hartree-Fock theory, or, den- 
sity functional theory with an approximate functional, 
to provide an accurate description. Exact diagonaliza- 
tion is prohibitively expensive because of the importance 
of Landau level mixing. Exact diagonalization wave func- 
tions within the LLL approximation are not sufficiently 
accurate on their own, but when they are multiplied by 
an optimized Jastrow and used as trial wave functions 




FIG. 5: Evolution of the ground state densities as a functions 
of magnetic field is increased for N = 6, obtained by VMC 
method, (a) Without Zeeman effect (b) Including Zeeman 
effect. 



in DMC, they lead to very accurate energies. The high 
accuracy of the fixed phase DMC energies shows that the 
phase of exact wave functions is not strongly affected by 
LL mixing. The method described is very efficient since 
the trial wave functions have determinantal coefficients 
that are independent of system parameters m* , ujq, B 
and e, and the Jastrow optimization process is very fast 
and stable. We have found that near the MDD-LDD 
transition, the pair densities for N — 4 have a well de- 
fined square structure, whereas those for N — 6 have 
little structure. When the Zeeman energy is included, 
the next ground state beyond the MDD state is fully po- 
larized and has L = N(N + l)/2. This leads to an abrupt 
redistribution of charge, as observed experimentally |(|. 
If instead the Zeeman energy is zero then several un- 
polarized and partially polarized states become ground 
states over a rather small range of magnetic fields and 
the charge redistribution is smaller and smoother. This 
could be verified experimentally. For instance, if a mate- 
rial for which the g— factor can be tuned to a small value 
is used, then the transition from unpolarized or partially 
polarized LDD states to fully polarized ones could be in- 
duced by increa sing the effective g— factor by using tilted 
magnetic fields |37| . 

The number of electrons that can be studied with this 
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method is limited by the number of Slater determinants 
in the LLL approximation. This increases rapidly as the 
angular momentum increases beyond the MDD value, 
L = N(N — l)/2). If the system is fully spin polarized, 
the number of determinants in the LLL approximation 
is much reduced and N can be much larger. When L is 
sufficiently large that the number of determinants in even 
the LLL becomes too large it is possible to dramatically 
reduce the number of determinants by employing com- 
posite fermion wave functions |3S|. which are very good 
approximations to the exact LLL wave functions. Once 



again Landau level mixing can be incorporated very effi- 
ciently with QMC. This has been done for a 15-electron 
dot in Ref . ElS 
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